% For "The Climate Risk Premium" (Lemoine, JAERE, 2020)

% sub_damagenodes_damexp:
% Hard-code quadrature nodes for the damage distribution, as calculated from Fortran90 version of truncated_normal_rule.

if abs(Params.mean_logdamage + 8.5564) < 0.1 % 8 nodes
    nodes_damage(:,1) =   [    -15.07278940310906
        -13.17492037994422
        -11.54883964855559
        -10.04903643527189
        -8.623128132509287
        -7.256247884446813
        -5.982134126277506
        -4.974056572177592      ];
    weights_damage(:,1) = [       0.3111344836753689E-04
        0.3340195662469874E-02
        0.5301025853248677E-01
        0.2371918858734920
        0.3806237918057968
        0.2456272758895616
        0.6988263839630796E-01
        0.1029284039151786E-01];
elseif abs(Params.mean_logdamage + 7.6771 ) < 0.1 % 8 nodes
    nodes_damage(:,1) =   [    -14.38357427482333
        -12.51836213934920
        -10.93058441120087
        -9.479722175512842
        -8.120464529897500
        -6.849863345844519
        -5.719372541659109
        -4.896434932418704     ];
    weights_damage(:,1) = [        0.1739687028997504E-04
        0.2040892314529469E-02
        0.3590114234047578E-01
        0.1828945421120232
        0.3491156959808180
        0.2876371609984619
        0.1161736756661191
        0.2621949371728249E-01 ];
elseif abs(Params.mean_logdamage + 6.5402 ) < 0.1 % 8 nodes
    nodes_damage(:,1) =   [     -13.55010013756096
        -11.73636844808782
        -10.20815093251008
        -8.831734709524101
        -7.570136387843832
        -6.431427672667297
        -5.474597163105261
        -4.833450050241599   ];
    weights_damage(:,1) = [      0.7051491837900867E-05
        0.9422890616186657E-03
        0.1923773849369622E-01
        0.1177517076653608
        0.2847243030041655
        0.3203301612085119
        0.1923468478455781
        0.6465990122923072E-01 ];
else
    error('Not hard-coded')
end

clear temp_value;